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ABSTRACT 

We report new magnetic field measurements for 14 classical T Tauri stars (CTTSs). We 
combine these data with one previous field determination in order to compare our observed field 
strengths with the held strengths predicted by magnetospheric accretion models. We use litera- 
ture data on the stellar mass, radius, rotation period, and disk accretion rate to predict the field 
strength that should be present on each of our stars according to these magnetospheric accretion 
models. We show that our measured field values do not correlate with the field strengths pre- 
dicted by simple magnetospheric accretion theory. We also use our field strength measurements 
and literature X-ray luminosity data to test a recent relationship expressing X-ray luminosity as 
a function of surface magnetic flux derived from various solar feature and main sequence star 
measurements. We find that the T Tauri stars we have observed have weaker than expected X-ray 
emission by over an order of magnitude on average using this relationship. We suggest the cause 
for this is actually a result of the very strong fields on these stars which decreases the efficiency 
with which gas motions in the photosphere can tangle magnetic flux tubes in the corona. 



Subject headings: accretion, accretion disks — line: profiles 
— stars: magnetic fields — stars: pre-main-sequence — 



stars: atmospheres — stars: formation 



1. Introduction 

It is now generally accepted that accretion of 
circumstellar disk material onto the surface of a 
classical T Tauri star (CTTS) is controlled by a 
strong stellar magnetic field (e.g. see review by 
Bouvier et al. 2006). It is within the disks around 
these low mass pre-main sequence stars that solar 
systems similar to our own form. Understanding 
the processes through which young stars interact 
with and eventually disperse their disks is critical 
for understanding the rotational evolution of stars 
and the formation of planets. A key question is 
to understand how young stars can accrete large 
amounts of disk material with high specific angu- 
lar momentum, yet maintain rotation rates that 
are observed to be relatively slow (e.g. Hartmann 
& Stauffer 1989, Edwards et al. 1994). This prob- 
lem is solved in current magnetospheric accretion 
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models by having the stellar magnetic field trun- 
cate the inner disk, typically near the corotation 
radius, and channel the disk material onto the stel- 
lar surface, most often at high stellar latitude. 

Support for magnetospheric accretion in CTTSs 
is significant. Current models can account for the 
relatively slow rotation of most CTTSs (Camen- 
zind 1990; Konigl 1991; Collier Cameron & Camp- 
bell 1993; Shu et al. 1994; Paatz & Camenzind 
1996; Long, Romanova, & Lovelace 2005). Studies 
of the spectroscopic and photometric variability of 
CTTSs have been interpreted in terms of magnet- 
ically controlled accretion (e.g. Bertout, Basri, & 
Bouvier 1988; Bertout et al. 1996; Herbst, Bailer- 
Jones, & Mundt 2001; Alencar, Johns-Krull, & 
Basri 2001) with the magnetic axis inclined to the 
rotation axis in some cases (e.g. Kenyon et al. 
1994, Johns & Basri 1995, Bouvier et al. 1999, 
Romanova et al. 2004). Models of high resolu- 
tion Balmer line profiles, computed in the context 
of magnetospheric accretion, reproduce many as- 
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pects of observed line profiles (Muzerolle, Calvet, 
& Hartmann 1998, Muzerolle et al. 2000). Such 
models can also reproduce aspects of the observed 
Ca II infrared triplet lines (Azevedo et al. 2006). 
T Tauri stars (TTSs) are observed to be strong 
X-ray sources indicating the presence of strong 
magnetic fields on their surfaces (see Feigelson & 
Montmerle 1999 for a review), and several TTSs 
have been observed to have strong surface mag- 
netic fields (Basri, Marcy, & Valenti 1992; Guen- 
ther et al. 1999; Johns-Krull, Valenti, & Koresko 
1999b, hereafter Paper I; Johns-Krull, Valenti, & 
Saar 2004, hereafter Paper II; Yang, Johns-Krull, 
& Valenti 2005, hereafter Paper III), and strong 
magnetic fields have been observed in the forma- 
tion region of the He I emission line at 5876 A 
(Johns-Krull et al. 1999a; Valenti & Johns-Krull 
2004; Symington et al. 2005), which is believed 
to be produced in a shock near the stellar surface 
as the disk material impacts the star (Beristain, 
Edwards, & Kwan 2001). 

Despite these successes, open issues remain. 
Most current theoretical models assume the stel- 
lar field is a magnetic dipole with the magnetic 
axis aligned with the rotation axis. However, re- 
cent spectropolarimetric measurements show that 
the fields on TTSs are probably not dipolar, either 
aligned with the rotation axis or not (Johns-Krull 
et al. 1999a; Valenti & Johns-Krull 2004; Daou, 
Johns-Krull, & Valenti 2005; Yang, Johns-Krull, 
& Valenti 2006). On the other hand, it is expected 
that even for complex magnetic geometries, the 
dipole component of the field should dominate at 
distance from the star where the interaction with 
the disk is taking place, so this may not contra- 
dict current theory. In the case of the Sun, the 
dipole component appears to become dominant 
at 2.5i?o or closer (e.g. Luhmann et al. 1998). 
For expected disk truncation radii of 3 - 10 i?* 
in CTTSs (see below and Table [1]), this suggests 
the dipole component will govern the stellar inter- 
action with the disk. Additionally, Gregory et al. 
(2006) show that accretion can occur from a trun- 
cated disk even when the stellar field geometry is 
quite complex; however, no study has considered 
the torque balance between a star and its disk in 
the case of a complex stellar field geometry. 

Magnetospheric accretion models for CTTSs 
have been developed by a number of investigators. 
Konigl (1991) first applied the work of Ghosh and 



Lamb (1979) to CTTSs, showing that an equi- 
librium state could exist if the stellar field was 
strong enough. The field would truncate the inner 
disk and the star would spin at the same angular 
velocity as the disk where it was truncated just 
inside the corotation radius. Additional authors 
have have studied this equilibrium state, utilizing 
different assumptions regarding the details of how 
the stellar field couples to the disk. Paper I exam- 
ined the models of Konigl (1991), Collier Cameron 
and Campbell (1993), and Shu et al. (1994) and 
presented equations from each work which predict 
the stellar magnetic field strength given the stel- 
lar mass (M„), radius (R#), rotation period (P ro t), 
and mass accretion rate (M). These three stud- 
ies each assume the stellar field is a dipole aligned 
with the rotation axis. In this simple case, there 
should be no photometric or spectroscopic vari- 
ability associated with accretion assuming the disk 
is uniform. As mentioned above, the observations 
clearly show substantial accretion variability, and 
they also show that the fields are probably not 
aligned dipoles, at least at the stellar surface. Nev- 
ertheless, it may well be that equilibrium models 
mentioned above give an indication of the mean, 
or time averaged behavior of a CTTS. 

Several authors have performed numerical sim- 
ulations of magnetospheric accretion. Typically, 
aligned stellar dipole fields are assumed and the 
simulations examine the instabilities that can re- 
sult between the star, its magnetosphere, and the 
disk. Many studies find that winds and/or jets are 
produced (e.g. Hayashi, Shibata, & Matsumoto 
1996; Goodson, Winglee, & Bohm 1997; Miller 
& Stone 1997; Goodson, Bohm, & Winglee 1999; 
Ferreira, Pelletier & Appl 2000; Romanova et al. 
2005). These studies usually examine only a few 
select cases and generally do not address how the 
rotation of the star may evolve: it is a fixed pa- 
rameter of the model. These simulations may be 
testable through variability studies (see Goodson 
et al. 1999). Recently, Long et al. (2005) present 
numerical simulations, again assuming an aligned 
dipole field geometry, with the aim of determin- 
ing if a time averaged equilibrium rotation rate 
is established. Despite variability in the accretion 
and wind flows, these authors find that the star 
attains an equilibrium rotation rate with a period 
equal to that of the disk where it is truncated in- 
side, but close to, the corotation radius. Long et 
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al. (2005) use the results of their simulations to 
derive a relationship for the equilibrium rotation 
period of CTTSs (their equation 12). This can 
be solved for the stellar magnetic field strength 
and results in an equation with the same depen- 
dence on the stellar (M*, P ro t) and accretion 
parameters (M) as found for Konigl (1991) and 
Shu et al. (1994) as given in Paper I. It should 
be noted however, that Matt and Pudritz (2004, 
2005a) have recently asserted that the spin down 
torque produced on the star in magnetospheric ac- 
cretion models may be up to an order of magni- 
tude lower than most previous investigators have 
calculated. Matt and Pudritz (2005) suggest that 
an accretion powered stellar wind is responsible 
for spinning down CTTSs. 

Given the theoretical uncertainty involved with 
magnetospheric accretion, we can try to use obser- 
vations to test whether the disk locking scenario 
proposed by the equilibrium models actually holds 
in CTTSs. One of the first attempts at this is the 
work of Stassun et al. (1999) who find no correla- 
tion between rotation period and the presence of 
an infrared (IR) excess indicative of a circumstel- 
lar disk in a sample of 254 stars in Orion. How- 
ever, IR excess alone is not a good measure of 
the accretion rate. Muzerolle, Calvet, and Hart- 
mann (2001) note that current theory predicts a 
correlation between rotation period and mass ac- 
cretion rate which they do not observe. Muzerolle 
et al. (2001) suggest that variations in the stellar 
magnetic field strength from star to star may ac- 
count for the lack of correlation. Indeed, Paper I 
emphasizes that there are several stellar and ac- 
cretion parameters that enter into the equilibrium 
relationship, and the stellar magnetic field remains 
the quantity measured for the fewest number of 
CTTSs. Due primarily to this lack of data, Johns- 
Krull and Gafford (2002) made the assumption 
that the field in fact does not vary significantly 
from star to star and they looked for correlations 
among the remaining stellar and accretion param- 
eters. The predictions based on aligned dipolc field 
geometries are not present in the data; however, 
if the dipole assumption is dropped, Johns Krull 
and Gafford showed that the predicted correla- 
tions are present in the samples they examined 
from CTTSs in Taurus. 

The success of the Johns-Krull and Gafford 
(2002) study depends on the constancy of the mag- 



netic field from one TTS to the next. From a dy- 
namo perspective, this may be a good assumption. 
For cool main sequence stars there is a well de- 
fined positive correlation between magnetic activ- 
ity indices (e.g. Ca II H&K emission, X-ray emis- 
sion) and inverse Rossby number (the convective 
turnover time divided by the rotation period) with 
saturation setting at large inverse Rossby num- 
ber (e.g. Vilhu 1984, Noyes et al. 1984). Johns 
Krull, Valcnti, and Linsky (2000) showed that due 
to their long convective turnover times and rela- 
tively short rotation periods, TTSs generally all 
lie in the saturated portion of this relationship. 
On the other hand, most TTSs are fully convec- 
tive or nearly so, so a solar-like interface dynamo 
(Durney & Latour 1978) is likely very inefficient 
or non-existent in these stars. Another possibility 
is a convective/turbulent dynamo such as origi- 
nally proposed by Durney, DcYoung, & Roxburgh 
(1993). Recent studies by Chabrier and Kuker 
(2006) and Dobler, Stix, and Brandenburg (2006) 
do predict that dynamo action in fully convective 
stars does correlate with rotation rate. A final 
potential contributor to the fields of TTSs is pri- 
mordial fields entrained in the star during the star 
formation process (e.g. Tayler 1987, Moss 2003). 
However, little is known about how these fields 
might vary from star to star in a given star forma- 
tion region. As a result, additional measurements 
of the magnetic field strength on a sizeable sample 
of CTTSs are required to test current magneto- 
spheric accretion models. 

To that end, we have been conducting obser- 
vations of TTSs at the NASA Infrared Telescope 
Facility (IRTF) to measure the Zeeman broaden- 
ing of several K band Ti I lines which are promi- 
nent in the spectra of cool stars. The first use of 
this line to measure a field on a TTS is in Pa- 
per I, where we lay out most of the basics of our 
observation and analysis technique. Details have 
been expanded upon somewhat in Papers II and 
III. In particular, in Paper III, it was shown that 
errors in effective temperature of ~ 200 K as well 
as errors of 0.5 dex in logg result in errors of less 
than 10% in the derived mean magnetic field on 
TTSs. This is because we are measuring actual 
Zeeman broadening of the photospheric absorp- 
tion lines instead of some secondary effect such as 
a change in line equivalent width. As a result, it is 
no longer necessary to perform the detailed anal- 
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ysis of optical high resolution spectra to obtain 
the best estimates of T c g and log g in order to ob- 
tain a good magnetic held measurement. What is 
needed is a good estimate of usini, and all the 
stars in our sample have reliable measurements 
from the literature. Here we focus on CTTSs for 
which we have the ancillary data required to pre- 
dict the stellar magnetic held strength from mag- 
netospheric accretion models. To do that, we use 
the equations (1-3) given in Paper I for the stud- 
ies of Konigl (1991), Collier Cameron and Camp- 
bell (1993), and Shu et al. (1994), and we use 
equation (12) of Long et al. (2005). In order to 
make the held predictions, we need stellar and ac- 
cretion parameters. Rotation periods come from 
Bouvier et al. (1993, 1995), with the exception 
of TW Hya which comes from Mekkaden (1998). 
Spectral type, stellar radius, and the mass accre- 
tion rate for most stars come from Gullbring et 
al. (1998). For DG Tau, DH Tau, and T Tau we 
take these values from White and Ghez (2001). In 
some cases (e.g. DH Tau), the spectral type is 
listed as spanning two subtypes (e.g. Ml-2) and 
here we average the subtypes (e.g. Ml. 5). For 
all stars, we use pre-main sequence evolutionary 
tracks of Siess, Dufour, and Forestini (2000) to de- 
rive the stellar mass. Our sample of stars is given 
in Table [T] Included in Table Q] is the observed 
X-ray luminosity for each star. These come from 
Neuhauser et al. (1995) except for DF Tau, DH 
Tau, and GM Aur which come from Damiani et 
al. (1995); and for TW Hya which comes from 
Kastner et al. (2002). The goal of this paper is to 
measure the mean magnetic held strength for all 
the stars in TableQ]so that the observations can be 
compared with the theoretical predictions. Some 
of the stars in our sample are known close bina- 
ries; however, in these cases one member of the 
binary dominates the optical and IR light, and it 
is this member to which the various stellar, accre- 
tion, and magnetic parameters apply. Included in 
our sample is TW Hya which we do not analyze 
here, but was analyzed in Paper III. In §2 we de- 
scribe our observations and data reduction. The 
magnetic held analysis is described in §3. In §4 we 
give our results and we compare our measured and 
the predicted field strengths. We also compare our 
measured field strengths and the observed X-ray 
luminosity of our targets to theoretical predictions 
based on main sequence stars. A discussion of our 



results is given in §5. 

2. Observations and Data Reduction 

All spectra were obtained with the CSHELL 
spectrometer (Tokunaga et al. 1990, Greene et al. 
1993) on the 3.0 m telescope at the NASA IRTF 
on Mauna Kea in Hawaii. Observations occurred 
during four observing runs: I) 10-12 December 
1996, II) 14-20 December 1997, III) 7-12 January 
2000, and IV) 23-28 November 2004 (see Table 
[2]). On the first run a l'.'O slit gave a FWHM of 
~ 4.7 pixels on the 256 x 256 InSb array detec- 
tor, corresponding to a spectral resolving power of 
R = A/AA ~ 21,500. For the second and fourth 
runs, a 0'.'5 slit yielded a FWHM of - 2.8 pixels 
on the detector giving R ~ 36, 000. For the third 
observing run, the same 0'.'5 slit yielded a spec- 
tral resolving power of only R ~ 27, 000. [Soon 
after run III, CSHELL was serviced to restore the 
spectral resolution.] CSHELL uses a continuously 
variable filter (CVF) to isolate individual orders 
of the echelle grating. Each star was observed in 
three settings, each covering ~ 0.0057 /xm. The 
first setting (1) contains two strong Ti I lines at 
2.2211 and 2.2233 //m. The second (2) contains 
two more strong Ti I lines at 2.2274 and 2.2311 
/im. The third setting (3) centered at 2.3120 /xm 
contains ~ 9 CO lines from the v = 2 — first 
overtone band. 

Each star was observed at two positions along 
the slit separated by 10". Multiple pairs of ex- 
posures were made of each CTTS. Total exposure 
time varied some from object to object and wave- 
length setting to wavelength setting. The mini- 
mum exposure time was 20 minutes and the max- 
imum was 1.5 hours. A typical total exposure time 
was 1 hour. Taking the difference of each image 
pair with the star moved along the slit removed 
detector bias, dark current, and the average back- 
ground due to night sky emission. Difference im- 
ages were flat fielded using a normalized spectrum 
of an internal continuum lamp that fully illumi- 
nated the slit. Stellar spectra were extracted fol- 
lowing the procedure described in Paper I, which 
includes the determination of an oversamplcd slit 
function and optimal extraction of the spectrum. 
Wavelengths were determined by htting nX as a 
function of line position for 7-10 lamp emission 
lines and then dividing by the order number, n, 
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of the desired spectral region. Calibration lines 
from several orders were obtained by changing the 
order-sorting CVF while keeping the grating fixed. 
In addition to CTTSs, each night we observed 
hot, rapidly rotating stars spanning a range of air- 
masses. Using these data, all spectra shown in this 
paper have been corrected for telluric absorption. 
For each CTTS in our sample, Table [2] gives a log 
of when the star was observed and which wave- 
length settings were covered. 

3. Analysis 

The most successful approach for measuring 
fields on late-type stars in general has been to 
measure Zeeman broadening of spectral lines in 
unpolarizcd light (e.g., Robinson 1980; Saar 1988; 
Valenti, Marcy, & Basri 1995; Johns-Krull & 
Valenti 1996; Papers I, II, & III). In the presence 
of a magnetic field, a given spectral line will split 
up into a number of components depending the 
atomic structure of the levels contributing to the 
line. Taking into account rotational and turbulent 
broadening, a strong magnetic field typically pro- 
duces a change in the shape of magnetically sen- 
sitive line profiles. There are several Ti I lines in 
the K band which are excellent probes of magnetic 
fields in late- type stars (e.g. Saar & Linsky 1985), 
and here we observe 4 of them (see §2). In addition 
to the 4 strongly Zeeman sensitive Ti I lines, our 
wavelength settings also cover a relatively weak, 
but detectable, Sc I line at 2.2265 /im which has 
reduced Zeeman sensitivity due to lower Lande- 
g values. We also observe several magnetically 
insensitive CO lines near 2.313 /im which serve 
as a check on our values of v sin i and turbulent 
broadening terms. Since we are observing CTTSs, 
there is some concern that these lines could form 
in the disk and trace disk kinematics instead of 
stellar properties. This issue was studied for CO 
by Casali and Eiroa (1996) and Johns-Krull and 
Valenti (2001). Both studies find that the K band 
CO lines of the vast majority of CTTSs arise in 
the stellar photosphere, and the CO lines observed 
below are fully consistent with formation on the 
star. Based on their behavior as a function of ef- 
fective temperature, the Ti I lines in the K band 
appear to trace similar temperature material as 
that traced by the CO (Wallace & Hinkle 1997). 
In addition, the limited studies of TTS which show 
no K band excess (Johns-Krull et al. 2004, Yang 



et al. 2005) reveal Ti I and CO line profiles with 
the same qualitative behavior as that found below: 
CO lines well matched by rotationally broadened 
stellar photosphere models and excess broadening 
in the magnetically sensitive Ti I lines. We there- 
fore assume that the line profiles are dominated 
by the stellar photosphere. 

The goal in this paper is to measure the mag- 
netic field on a sample of CTTSs by modeling the 
profiles of their K band photospheric absorption 
lines. Our spectrum synthesis code and detailed 
analysis technique for measuring magnetic fields 
on TTSs is described in Papers I - III. Here, we 
simply review some of the specific details rele- 
vant to the results shown here. In order to syn- 
thesize the stellar spectrum, we must first spec- 
ify the atmospheric parameters: effective temper- 
ature (Toff), gravity (log g), metallicity ([M/H]), 
microturbulence (v m i c ), macroturbulence (w mac ), 
and rotational velocity (vsini). Following Papers 
I — III, a fixed value of 2 km s _1 is adopted 
for w mac . Valenti, Piskunov, & Johns-Krull 1998 
found that microturbulence and macroturbulence 
were degenerate in M dwarfs, even with very high 
quality spectra. Therefore, microturbulence is 
neglected here, allowing v mac to be a proxy for 
all turbulent broadening. Solar metallicity is as- 
sumed for all stars. We use v sin i values measured 
using the CO spectra presented here by Johns- 
Krull and Valenti (2001) for all the target stars an- 
alyzed in that paper. We adopt vsini values from 
Hartmann et al. (1986) for CY Tau and CM Aur, 
and we adopt the v sin i value of Basri and Batalha 
(1990) for DG Tau. The T cS is based on spectral 
type using the calibration of Johnson (1966) and 
log g is estimated by placing objects in the HR 
diagram. This typically produces uncertainties of 
< 200 K in T e ff and < 0.5 dex in log g which trans- 
late into uncertaities of 10% or less in the derived 
magnetic field values (Paper III). Once we have 
estimates of T e ff and log g, we use the "next gen- 
eration" (NextGen) model atmosphere (Allard & 
Hauschildt 1995) from the published grid which is 
closest to these values. In the range of interest 
here, the NextGen models are on a grid every 100 
K (e.g. 3900, 4000, 4100) in T cff and every 0.5 dex 
(e.g. 3.5, 4.0, 4.5) in log g. For BP Tau, we use the 
NextGen model with T c g= 4100 K which is closest 
in temperature to the value we derived in Paper I 
(its spectral type of K7 suggests T e s= 4000 K). 
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Lastly, we solve for the K band veiling in each 
star along with the stellar magnetic field proper- 
ties. We initially decided to fit for a single veiling 
parameter which would apply to all three spectral 
regions we observe in each star. This turned out 
to not be practical, as a few stars could not be well 
fit by a single veiling value. The different wave- 
length regions for our stars were observed with 
a time interval between settings of several hours 
up to a couple of days. Classical TTSs regularly 
show significant variations in their K band flux 
on timescales as short as a day (and occasionally 
shorter), likely as a result of accretion variabil- 
ity (Skrutskie et al. 1996, Eiroa et al. 2002). We 
therefore decided to let the veiling be a free param- 
eter for each star in each of the three wavelength 
settings analyzed here. 

The free parameters of our fit to the observed 
spectra are then the value of the K band veiling 
in each wavelength region (this is the only free pa- 
rameters in the CO wavelength region) and the 
magnetic field properties of the star. It was found 
in Papers I - III that the spectra of the TTSs stud- 
ied there could not be fit with a single value of the 
magnetic field strength. Instead, a distribution of 
magnetic field strengths is required. It was also 
found that fits to the spectra are degenerate in 
the derived field values unless we limit the fit to 
specific values of the magnetic field strength, sep- 
arated by about 2 kG, which is the approximate 
"magnetic resolution" of our data. Therefore, we 
use the same limitations when fitting the data pre- 
sented here. We assume the star is composed of 
regions of 0, 2, 4, and 6 kG magnetic field, and 
we solve for the filling factor of each of these com- 
ponents. The different regions are assumed to be 
well mixed over the surface of the star - different 
components are not divided up into well defined 
spots or other surface features. Another key as- 
sumption (which will be discussed further in §5) is 
that the temperature structure in all the field re- 
gions is assumed to be identical for the purpose of 
spectrum synthesis: again, the fields are not con- 
fined to cool spots or hot plage-like regions. In 
total then, we fit for 6 free parameters: the 3 veil- 
ing values and the filling factor of the 3 non-zero 
field regions (the sum of the filling factors must 
be 1.0). In Figures [T] and [2] we show our spectral 
fits to two of our CTTSs, the first (DK Tau) with 
a high average surface field, and the second (DE 



Tau) with a low average surface field. We char- 
acterize the field on all our CTTSs by computing 
the mean magnetic field from our fits: B = S-Bj/j 
where Bi is the value of the field (0, 2, 4, 6 kG) 
in each fitted region and fa is the filling factor of 
these field components. The formal uncertainty in 
the fit to all our stars is quite small. The true un- 
certainty is dominated by uncertainties in the field 
resulting from errors in our adopted temperature 
and gravity, which as discussed above are ~ 10%. 
We therefore adopt this value for the uncertainty 
in our mean magnetic field for each star. 

While DK Tau's Ti I lines show obvious broad- 
ening relative to the nonmagnetic model in Figure 
[TJ the situation is a little less clear for DE Tau in 
Figure [2] Blowing up the plot, it is apparent that 
all the Ti I lines are broader than the null field 
model which is a good match to the width of the 
CO lines. This extra broadening in the Ti I lines 
can be quantified by a cross correlation analysis 
as well. We construct a kernel which is the zero 
field model but with a wsinz = 1.0 km s _1 in- 
stead of 7 km s _1 which is the measured value for 
DE Tau. We then cross-correlate this kernel with 
the observed DE Tau spectrum and fit the result 
with a Gaussian profile to determine the FWHM 
of the cross correlation function. Doing this on 
the CO portion of the spectrum (lower panel of 
Figure [2]), we find a FWHM= 6.35 pixels. Cross 
correlating the 1.0 km s _1 kernel with the actual 
fit to the CO lines (dashed line of Figure [5] where 
vsini= 7.0 km s" 1 ), we find the FWHM= 6.54 
pixels (which is slightly greater - by 0.19 pixels 
- than for the observed spectrum). We can esti- 
mate the range of FWHM we should measure by 
taking the fit to the spectrum shown in the dashed 
line and adding random noise to it to simulate the 
signal-to- noise ratio (S/N) of our observed spec- 
trum (115 for DE Tau). We then cross correlate 
this synthetic data with the same kernel above and 
repeat this experiment 100 times. Doing so, we re- 
cover the a mean FWHM of 6.54 pixels and find 
a standard deviation in our recovered FWHM val- 
ues of 0.09 pixels. Thus, the FWHM we recover 
from the observations of the magnetically insensi- 
tive CO lines is only 2.1<r less than that recovered 
from the CO lines of the null field model for DE 
Tau. When we perform the same analysis on the 
magnetically sensitive Ti I lines, we find that the 
FWHM of the cross correlation function using the 
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observations is 9.97 pixels, whereas the null field 
model gives a FWHM= 8.66 pixels with a stan- 
dard deviation of 0.13 pixels. Thus the observed 
magnetically sensitive Ti I lines give a width lO.ler 
broader than does the null field model. Using the 
actual magnetic fit to DE Tau and again adding 
100 different realizations of random noise to it to 
create synthetic data for the cross correlation anal- 
ysis yields a FWHM= 10.12 ± 0.17 pixels, which 
is within 0.9<r of that recovered from the observed 
spectrum. We are therefore confident in the excess 
width of the Ti I lines in DE Tau (B = 1.12 kG) 
and in the other "low" field measurement stars 
(see below) . The analysis of Paper III shows that 
even for these low mean field values, the uncer- 
tainty due to temperature and gravity errors again 
lead to a mean field uncertainty of ~ 10%, which 
is the field uncertainty adopted here. 

4. Results 

The results of our spectral fitting are given in 
Table [3l Reported here are the derived K band 
veiling, r^, the surface averaged magnetic field, 
B, for each star, and a prediction for the X-ray 
luminosity based on the measured field values (dis- 
cussed below). One star in our sample, GM Aur, 
was observed in only one wavelength setting: the 
Ti I setting (2) described above. Results for this 
star are somewhat less firm that for the others. 

4.1. K Band Veiling 

We determine 3 estimates of the K band veil- 
ing (rx) for each star (one from each wavelength 
region) which we can compare against each other 
and against previous veiling estimates. In Table 
[3] we tabulate the mean value of rx determined 
from the 3 settings and give as an uncertainty es- 
timate the standard deviation of the 3 values. 
These veiling estimates generally agree with previ- 
ous studies (Folha & Emerson 1999, Johns-Krull 
& Valenti 2001). We note though that by solv- 
ing for the veiling in each wavelength region sepa- 
rately, our magnetic field results are free of system- 
atic effects caused by differences in line strength 
between the actual star and our spectrum synthe- 
sis. Such differences are adjusted for by changes 
in the derived r^- values. It is the shape of the 
Ti I lines that constrains the magnetic field vales 
we derive. In support of this idea, we looked for a 



correlation between the veiling values found here 
and the mean fields measured. We also looked for 
a correlation between the difference in the veil- 
ing found here compared to in Johns-Krull and 
Valenti (2001) and the mean field values. No cor- 
relation is found in either case. 

For one star, T Tau, our derived r^ is sub- 
stantially different than found in previous studies. 
At K0, T Tau is the earliest spectral type star in 
our study, and one of the earliest spectral type 
CTTS known. The veiling of Johns-Krull and 
Valenti (2001) is derived based on a later spectral 
type template star. Similarly, Folha and Emer- 
son (1999) report a veiling of rx — 2.5 based on 
comparison with a later type template. The Ti I 
and CO lines used here get much weaker in the 
observed spectra of dwarfs and giants as you go 
from K7 to K0 (e.g. Wallace & Hinkle 1997). This 
difference in intrinsic line strength will artificially 
enhance the derived veiling for a K0 star when us- 
ing a later template. The veiling we derive here 
is quite low, which might be the result of the line 
strength in our synthetic spectra being too low as 
a result of some error in the line data or the model 
atmosphere. To test this, we compared our model 
spectra for T Tau to observed spectra of the K0 
dwarf, HR 166. This star was observed in setting 
(2) during the 4th observing run (see §2). Valenti 
and Fischer (2005) report T cS = 5221 K and [M/H] 
= 0.16 for this K0V star. We observe the lines in 
HR 166 to be quite weak, and the 2.2311 (im Ti I 
line is strongly affected by telluric absorption, so 
we ignore it here. We measure an equivalent width 
of W eq = 49.4 ± 4.8 mA for the 2.2274 /im based 
on the 4 spectra observed of HR 166. Our kG 
synthetic spectrum for T Tau gives an equivalent 
width of W eq = 58.0 mA, in good agreement with 
our observed value. The fact that the synthetic 
spectrum has a larger equivalent width means that 
we may actually be over estimating the veiling, 
not underestimating it; however, the correction is 
small (~ 0.2). Therefore, we suggest our veiling 
measurement for T Tau is reliable, and that previ- 
ous estimates are too large as a result of the later 
type template stars used for comparison. 
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4.2. Magnetic Field Properties 

4-2.1. Comparison with Magneto spheric Accre- 
tion Theory 

Tabic [3] gives the mean magnetic field values, 
B, for each star analyzed here. These can be com- 
pared directly with the predicted field strengths 
given in Table [TJ Figure [3] compares the values of 
B we measure here with the predictions of magne- 
tospheric accretion theory. As mentioned above, 
all 4 field predictions given in Table [1] are very well 
correlated with each other. From the standpoint of 
looking for a correlation between the observed and 
predicted fields, it does not matter which set of 
theoretical predictions we use. Figure [3] plots the 
predictions based on the work of Shu et al. (1994). 
As is quite evident from a visual inspection of the 
figure, there does not appear to be a correlation 
between the observed and predicted field values. 
We can quantify this by computing the linear cor- 
relation coefficient (also called Pearson's r value) 
and its significance (e.g. Press et al. 1986). Com- 
puting this for the data in Figure 02 we find a 
correlation coefficient of r = 0.08 which has an as- 
sociated false alarm probability of 79%, indicating 
no correlation in the data at all. While there is no 
apparent correlation in Figure [3j the good news is 
that the observed fields generally lie on or above 
the line indicating equality. Thus, the measured 
fields on these CTTSs are approximately the right 
magnitude for magnetospheric accretion to work. 
We will return to this point further in §5 below. 

When examining Tables Q] and [3] and Figure 
[3j worry may be aroused by the specific choice 
of mass accretion rates used. It is well known 
that mass accretion rates for CTTSs are difficult 
to measure, and estimates for specific stars can 
vary by an order of magnitude or more (see dis- 
cussion in Gullbring et al. 1998). Fortunately, the 
predicted field strengths vary only as the square 
root of the mass accretion rate ( Johns-Krull et al. 
1999b), making them somewhat less sensitive to 
the problems associated with measuring the accre- 
tion rate. In addition, different studies generally 
do find good correlation of their measured accre- 
tion rates from star to star, even if they differ on 
the overall magnitude (again, see Gullbring et al. 
1998). Therefore, we have some confidence that 
the variation of the predicted fields from star to 
star can be reasonably well estimated. In an ef- 



fort to see how the results presented here depend 
on the stellar and accretion parameters used for 
the stars studied, Table EH gives the mass accre- 
tion rates and predicted fields for 5 relatively large 
studies which contain many of the stars observed 
here. Since the derived mass accretion rates de- 
pend to some extent on the mass and radius as- 
signed to the star, in Table |4] the predicted fields 
are calculated using the mass and radius adopted 
for each star in the relevant study. The studies 
presented in the table are Valenti et al. (1993, la- 
belled VBJ), Hartigan et al. (1995, labelled HEG), 
Gullbring et al. (1998, labelled GHBC), Calvet 
and Gullbring (1998, labelled CG), and White and 
Ghez (2001, labelled WG). While accretion rates 
can vary by over an order of magnitude bewteen 
these different studies, there is generally good cor- 
relation in the predicted fields based on the dif- 
ferent studies, particularly when allowing removal 
of a single obviously discrepant star. For exam- 
ple, comparing field predictions from HEG with 
those of GHBC produces a correlation (r = 0.73) 
with a false alarm probability of 1% even though 
the accretion rates differ by an order of magnitude 
in many cases. Comparing field predictions from 
VBJ with GHBC gives a correlation (r = 0.74) 
with a false alarm probability of 2% (due to the 
smaller number of stars in common) , and when re- 
moving the obviously discrepant GI Tau, the cor- 
relation improves to r = 0.99 with a false alarm 
probability of 2.7 x 10~ 6 . 

On the other hand, none of these studies pro- 
duced a good correlation with the observed mean 
field strengths. The bottom line of the table 
gives the correlation coefficient and associated 
false alarm probability when comparing the pre- 
dicted fields from each study with the observed 
fields. Plots of the observed versus predicted fields 
appear similar to Figure [3l showing a flat relation 
with no apparent correlation even when the most 
discrepant point is ignored. While, mass accre- 
tion rate is the parameter that varies the most for 
each star from one study to the next, the predicted 
fields turn out to be very sensitive to the adopted 
value of the stellar radius. The predicted field 
strengths vary as i?J (Johns-Krull et al. 1999b), 
which is usually determined by measuring the stel- 
lar luminosity and effective temperature. Effective 
temperatures usually come from spectral type and 
an adopted spectral type - effective temperature 
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calibration; however, this has a number of its own 
problems (e.g. Padgett 1996, Huerta 2007) and is 
deserving of more attention. 

4- 2. 2. Relationship to X-ray Emission 

Stellar magnetic fields are also believed to give 
rise to "activity" in late- type stars. Activity is 
typically traced by line emission or broad band 
emission at high energy wavelength such as X- 
rays. Pevtsov et al. (2003) study the relation- 
ship between magnetic flux and X-ray luminosity 
on the Sun and cool stars, generally finding excel- 
lent correlation between these two quantities rang- 
ing over almost 11 orders of magnitude in each 
quantity. Pevtsov et al. (2003) include prelim- 
inary measurements for 6 of the CTTSs shown 
here taken from Johns-Krull and Valenti (2000) 
and Johns-Krull et al. (2001). We use the X-ray 
luminosity-magnetic flux relationship of Pevtsov 
et al. (2003) with our magnetic field measure- 
ments and the stellar radii given in Table [T] to cal- 
culate the expected X-ray luminosity for each of 
our stars. These predicted X-ray values arc given 
in the last column of Table [3] Figure |4] plots the 
measured X-ray luminosities from Table [T] versus 
these predicted X-ray values. All but one of the 
CTTSs in our sample show lower X-ray emission 
than would be predicted by their magnetic field 
properties, and many of them are low by more 
than an order of magnitude. Only TW Hya shows 
a higher level of measured X-ray emission relative 
to the prediction, and perhaps significantly, this 
CTTSs is the star closest to the main sequence in 
our sample due to its age of ~ 10 Myr. We discuss 
this issue further in the next section. 

4-2.3. Field Measurements and Pressure Equipar- 
tition 

From the perspective of cool star research, the 
high apparent field strengths on these CTTSs is 
a little surprising. Spruit & Zweibel (1979) com- 
puted flux tube equilibrium models, showing that 
magnetic field strength is expected to scale with 
gas pressure in the surrounding non-magnetic pho- 
tosphere. Similar results were found by Rajaguru, 
Kurucz, and Hasan (2002) who find that convec- 
tive collapse of thin magnetic flux tubes will set 
a field strength limit of ~ 750 G in the photo- 
sphere of a star with T c g= 4000 K and log g= 3.5. 
Convective collapse is a mechanism for forming 



partially evacuated flux tubes, and Rajaguru et 
al. (2002) find evidence that this mechanism may 
not operate effectively at low T e g; however, they 
find a maximum field strength of ~ 1300 G if 
the field is to remain in pressure equilibrium with 
the surrounding unmagnetized photosphere. Field 
strengths set by such pressure equipartition con- 
siderations appear to be observed in G and K 
dwarfs (e.g. Saar 1990, 1994, 1996) and possibly 
in M dwarfs (e.g. Johns-Krull & Valenti 1996). 

TTSs have relatively low surface gravities and 
hence low photospheric gas pressures, so that 
equipartition flux tubes would have relatively low 
magnetic field strengths compared to cool dwarfs. 
Indeed, Safier (1999) examined in some detail the 
ability of TTS photospheres to confine magnetic 
flux tubes via pressure balance with the surround- 
ing quiet photosphere, concluding that the maxi- 
mum field strength allowable on TTSs is substan- 
tially below the few detections reported at that 
time. The maximum field strength allowed for a 
confined magnetic flux tube is B eq — (SnPg) 1 / 2 
where P g is the gas pressure at the observed level 
in the stellar atmosphere. Here, we take as a lower 
limit (upper limit in pressure) the level in the at- 
mosphere where the local temperature is equal to 
the effective temperature. This is the approximate 
level at which the continuum forms, with the Ti I 
lines forming over a range of atmospheric layers 
above this level. Our computed values of B eq are 
given in Table O In Figure [5] we plot our observed 
mean magnetic field strengths versus these B eq 
values for our entire sample. In all cases, the ob- 
served mean magnetic field is greater than the field 
predicted by pressure equipartition arguments. In 
many cases the observed mean field exceeds B eq 
by a factor of more than 2. 

4-2.4- Origins of the Strong Fields 

Lastly, we compare our mean field measure- 
ments with various stellar parameters in an effort 
to uncover the origin of the strong magnetic fields 
observed on these stars. We compare both the 
mean magnetic field strength, B, and the total 
magnetic flux, $, to a number of stellar parame- 
ters, including those thought to be important in 
the dynamo generation of magnetic fields. The re- 
sults of these comparisons are summarized in Ta- 
ble El Given in this table is the Pearson linear 
correlation coefficient, r, and its false alarm prob- 
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ability, f p , which measures the degree to which the 
quantities being compared are correlated with one 
another. As Table [5] shows, there is no significant 
correlation of B or $ with any stellar parameter 
or common dynamo variable, including the con- 
vective turnover time, r c [values of r c were kindly 
provided by Y.-C. Kim and were computed using 
models described in Kim and Demarque (1996) 
with input physics described in Preibisch et al. 
(2005)]. The most significant correlation found 
is <fr versus L*; however, we must be careful with 
such a correlation due to the common role the stel- 
lar radius plays in both variables. The magnetic 
flux is <f> = BRl, and the stellar radius is calcu- 
lated from the measured L* and T e ff. Thus, we 
expect to find a correlation, and when comparing 
i?* versus L* we find r — 0.73 and f p = 0.002. 
This correlation is stronger than that found us- 
ing $, suggesting the correlation of $ with is 
driven entirely by the stellar radius, with the mag- 
netic field playing no role. Interestingly, the next 
most significant correlation is $ versus r c ; how- 
ever, this correlation is again likely driven by the 
correlation of i?* with r c (r = — 0.63, f p = 0.01), 
with the magnetic field playing no role. Within 
this sample of 15 CTTSs, we then find no corre- 
lation of the measured magnetic field values with 
any stellar or dynamo parameters. 

5. Discussion 

We provide new measurements of the mean 
magnetic field strength on 14 CTTSs. The de- 
tected fields range from 1.12 to 2.90 kG. These 
measurements combined with previous reports of 
magnetic fields on TTSs (Basri et al. 1992; Guen- 
ther et al. 1999; Johns-Krull et al. 1999b, 2004; 
Yang et al. 2005) suggest that the strong majority 
of TTSs are covered by kilogauss magnetic fields. 
As found in Papers I - III, we use a distribution 
of magnetic field strengths to fit the observed K 
band Ti I profiles. We model this distribution as if 
it were a spatial distribution on the stellar surface: 
the model assumes distinct regions with different 
field strengths. A concern is that the distribution 
is in fact a distribution with depth in the magnetic 
regions. If this were the case, we would be overes- 
timating the total filling factor of magnetic regions 
as well as the mean field strength. Observations 
of sunspots in the K band show that the Ti I lines 
in this spectral region display profiles characteris- 



tic of a single magnetic field value (e.g. Wallace 
& Livingston 1992), despite the likely variation of 
the field with depth in a sunspot (e.g. Mathew et 
al. 2003). Therefore, we believe the field distribu- 
tions found in Papers I - III and used here repre- 
sent actual spatial distributions of field strengths 
on the stellar surface. 

The observed mean field values for the CTTSs 
are all larger than the field strength for these stars 
predicted by pressure equilibrium arguments (Fig- 
ure [5]). This suggests that the magnetic pressure 
dominates the gas pressure, in some cases by a 
significant amount, in the photospheres of these 
stars. If this is the case, there should be no field 
free regions on these stars, though we do use a 
kG field component when fitting their spectra. 
This issue is discussed in §4.2.3 and in Papers I 
and II. In these papers, it was found that using fits 
with 0, 2, 4, and 6 kG components gave essentially 
identical mean field measurements as those from 
fits using 1, 3, 5, and 7 kG components. Thus, 
the data do not demand a field free region, and 
such a region may not exist on these stars. With 
strong magnetic field essentially covering the en- 
tire surface of TTSs, we expect to see some effect 
of the fields on the atmosphere, perhaps analo- 
gous to the fields of sunspots inhibiting convec- 
tion (e.g. Foukal 1990). Such an effect might be 
diagnosed through the study of line bisectors (e.g. 
Gray 1982); however, little work in this area has 
been done for late K spectral types, so the intrin- 
sic bisector shape for stars of this spectral type is 
unknown. 

On the other hand, we may have evidence for 
the feedback of these strong fields on the star 
through the X-ray activity displayed by TTSs. 
Feigelson et al. (2003), find that X-ray emission 
in TTSs saturates on average at a level approx- 
imately an order of magnitude below that dis- 
played by main sequence stars. In addition to 
this apparent global reduction in X-ray emission 
of TTSs, observations also reveal that accreting 
young stars as a group show X-ray emission levels 
that are reduced by a factor of 2-3 relative to non- 
accreting young stars, even after dependencies on 
mass and age are taken into account (Flaccomio 
et al. 2003, Stassun et al. 2004, Preibisch et al. 
2005, Feigelson et al. 2007). At least two expla- 
nations have been advanced that can explain this 
difference between stars with and without accre- 
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tion signatures. Preibisch et al. (2005) suggest 
that the higher density of mass loaded accreting 
field lines prevents these portions of CTTSs mag- 
netospheres from heating up to temperatures re- 
quired for X-ray emission. Jardine et al. (2006) 
present a model of TTS cornonae with realistic 
levels of complexity and show that this can pro- 
duce a similar range of X-ray emission levels as 
that which is observed in for example, the Chan- 
dra Orion Ultradeep Project (COUP). Jardine et 
al. (2006) suggest that the presence of a disk, par- 
ticularly in lower mass CTTSs, acts to strip the 
outer parts of the stellar corona, thereby reducing 
the total X-ray emission. This naturally predicts 
an anticorrelation between X-ray emission and the 
presence of a disk, at least in the lower mass young 
stars such as those studied here. While there is no 
observed correlation of X-ray emission with the 
presence or absence of an inner disk as diagnosed 
by K-band emission for the COUP dataset (e.g. 
Feigelson et al. 2007), this coronal stripping may 
be related to the action of accretion. While these 
ideas may explain the factor of 2-3 difference in 
accreting versus non-accreting young stars, they 
do not explain the order of magnitude reduction 
in X-ray emission for TTSs as a whole noted by 
Feigelson et al. (2003). Here, we find that the ob- 
served X-ray emission from our sample of CTTSs 
is substantially below what is expected given their 
strong magnetic field values (Figured]), often by 
more than an order of magnitude. While the ori- 
gin of coronal heating in the Sun and solar like 
stars is still debated, most theories assume it is due 
in part to convective motions moving and buffet- 
ing photospheric magnetic flux tubes (e.g. Fisher 
et al. 1998, Foukal 1990). We suggest that the 
strong magnetic fields which appear to cover the 
surface of TTSs reduce the efficiency with which 
convective motions can stress these stellar mag- 
netic fields, which in turn results in reduced coro- 
nal heating relative to what would otherwise be 
expected. 

The original motivation for this study was to 
test magnetospheric accretion models for CTTSs. 
Figure [3] shows the relationship between magnetic 
field strengths predicted from these models com- 
pared with the measured field strengths. Clearly, 
there is no correlation. Magnetospheric accretion 
models typically assume an aligned dipolar mag- 
netic field at the surface of a CTTS. It is now 



fairly clear that the fields of TTSs are not dipo- 
lar (Johns-Krull et al. 1999a, Valenti & Johns- 
Krull 2004, Daou et al. 2006, Yang et al. 2006). 
The magnetic geometry must be complicated at 
the surface of TTSs, as on the Sun and all other 
cool stars. However, the dipolar component may 
dominate at the inner edge of the accretion disk, 
several stellar radii from the surface. This notion 
is supported by observations of circular polariza- 
tion in the narrow component of the He I emission 
line which forms near the base of the accretion 
flows where the material strikes the stellar surface 
(Johns-Krull et al. 1999a). Time series polariza- 
tion in this line shows smooth variations charac- 
terized by few if any polarity reversals, and are 
consistent with the idea that the disk is interact- 
ing with a dipole-like field geometry (Valenti & 
Johns-Krull 2004, Symington et al. 2005, Yang et 
al. 2007). 

Magnetospheric accretion models use an as- 
sumed dipolar field geometry to predict which stel- 
lar magnetic field lines couple the stellar surface 
to the disk. Since the actual magnetic geometry is 
more complicated, the fraction of the stellar sur- 
face with field lines coupled to the disk (including 
those along which accretion will occur) can differ 
significantly from predicted values (see Gregory et 
al. 2006) for a discussion of this point. This means 
the relationships given in Paper I and used to cal- 
culate the predicted fields given in Table Q] are ac- 
tually underconstrained, even though we have now 
measured B at the surface. However, Johns-Krull 
and Gafford (2002) point out that a more general 
relationship for the Shu et al. (1994) model can 
be determined which does not rely on the assump- 
tion of a dipolar field geometry. This formulation 
characterizes the magnetic coupling in terms of the 
surface flux that threads the disk outside corota- 
tion, rather than the equatorial field strength of an 
assumed dipole. Valenti and Johns-Krull (2004) 
use the results of Johns-Krull and Gafford (2002) 
to derive an equation for the filling factor of accre- 
tion zones, f a , on CTTSs [their equation (1)]. We 
use our magnetic field results with this equation 
to predict f a for each of our stars. These accretion 
filling factors range from 0.57% to 3.09%, in good 
agreement with the range of observed f a values 
(Valenti, Johns, & Basri 1993; Calvet & Gullbring 
1998). The predicted f a values range over a fac- 
tor of 5.4 while the observed mean field strengths 
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range over a factor of 2.6. In the original mag- 
netospheric accretion models, f a is essentially as- 
sumed constant and the predicted field strengths 
range over a factor of 8 (Table[T]). We suggest that 
f a is more likely than B to vary by a large factor. 
Our measurements of B support this conclusion. 
Additionally, one could look for a correlation be- 
tween the fields predicted in Table [T] and the fields 
anchored in the accretion flows and diagnosed by 
the He I line. Johns-Krull et al. (2001) do this, 
finding a good correlation; however, we caution 
that they only had a sample of 4 stars. 

Given that the magnetic fields on the surface 
of TTSs may not be dipolar, one question that 
comes to mind is whether the observed fields on 
CTTSs are indeed strong enough to enforce disk 
locking in these young stars. The magnetosphcric 
accretion models require some (though admittedly 
uncertain) field strength in the disk, several stel- 
lar radii above the star, to enforce disk locking. 
The dipole assumption can then be used to predict 
what field strength is implied at the stellar surface. 
If the field geometry is much more complex than 
this, a larger field strength will be required since 
the higher order field components have a stronger 
dependence on radius than a dipole does. The con- 
tribution of the various field components at the in- 
ner disk edge is uncertain, but as mentioned above 
the He I observations suggest the dipole compo- 
nent may dominate. This field component is un- 
certain observationally, but it can be constrained 
through spectropolarimetry. Smirnov et al. (2003) 
report a marginal detection of a net field of 150 G 
on T Tau which was not confirmed by Smirnov 
et al. (2004) or Daou et al. (2006). Yang et al. 
(2007) detect a net field of 149 G on TW Hya on 
one night of their 6 night monitoring campaign 
on this star, finding only upper limits of ~ 100 
G on the other nights. Other studies also only 
find upper limits (3cr) of 100-200 G (Johns-Krull 
et al. 1999a, Valenti & Johns-Krull 2004). The 
relationship between the net field observed spec- 
tropolarimetrically and the dipole component of 
the star depends on knowing the angle between 
the line of sight and the dipole axis; however, the 
general lack of detection suggests the dipole com- 
ponent on TTSs is small. For example, a dipole 
with a 3 kG polar field strength (as suggested for 
BP Tau by Symington et al. 2005) observed at 
an angle of 45° would produce a net field of 690 



G, well above current detection limits. Accretion 
hot spots are expected to have essentially no effect 
on net field measurements from spectropolarime- 
try (Smirnov et al. 2005) since the acretion filling 
factors are typically ~ 1% (Valenti et al. 1993, 
Calvet & Gullbring 1998). 

Compared with at least some of the predicted 
values in Table [TJ the general lack of strong po- 
larimetric field detections in TTSs might suggest 
the fields are actually too weak to enforce disk 
locking in these stars. However, the general agree- 
ment in magnitude between the predicted and ob- 
served accretion filling factors described above and 
in Valenti and Johns-Krull (2004) suggests that 
the trapped flux model of Shu et al. (1994, see also 
Ostriker & Shu 1995) can still enforce disk locking 
with the field strengths observed here. Indeed, 
Johns-Krull and Gafford (2002) find that this 
trapped flux model predicts correlations among 
stellar and accretion parameters that are actually 
observed in the data, whereas correlations based 
on assuming a dipole field geometry are not ob- 
served in the data. As sectropolarimetric field 
measurements get more precise, and the dipole 
component gets more constrained, comparisons 
with all the field predictions in Table Q] will be- 
come more critical. In addition, more theoretical 
work is called for to explore whether realistic mag- 
netic field geometries and strengths can produce 
disk locking in magnetospheric accretion models. 
First steps in exploring the role of realistic field 
geometries in CTTSs accretion have been made 
by Gregory et al. (2006); however, torque balance 
calculations still need to be made for these types 
of models. 

Lastly, we consider the origin of the strong fields 
observed in TTSs. Chabrier and Kuker (2006) 
assert that primordial, or fossil, magnetic fields 
in fully convective stars can only survive for a 
timescale of ~ -R*/??, where is the stellar 
radius and rj is the turbulent magnetic diffusivity. 
The value of works out to be on the order of 10 
- 100 years. Thus, Charbier and Kuker (2006) as- 
sert that the fields of fully convective stars, such as 
most TTSs, must be produced by a turbulent dy- 
namo. On the other hand, Tayler (1987) and Moss 
(2003) suggest that primordial fields can survive 
much longer in pre-main sequence stars. Chabrier 
and Kuker (2006) and Dobler et al. (2006) both 
find that dynamo action in fully convective stars 
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should correlate with rotation rate. Motivated by 
this, we looked for correlations between the ob- 
served magnetic fields and various stellar param- 
eters including rotation (see §4.2.4). No correla- 
tions were found, and indeed stars like CY Tau and 
DH Tau show very similar stellar properties (T c g, 
mass, radius, rotation rate, convective turnover 
time) but have mean fields that differ by over a 
factor of 2.3. In addition, the work of Chabricr 
and Kiiker (2006) and Doblcr et al. (2006) gener- 
ally find magnetic field strengths somewhat lower 
than the equipartition values, while the observed 
fields are significantly stronger than this. Ad- 
ditional constraints on fully convective dynamo 
models can also be made by comparing predicted 
field geometries with those derived from Zeeman 
Doppler imaging of fully convective stars. Donati 
et al. (2006) do this for an M dwarf, finding incon- 
sistencies in the observed field geometry relative 
to the predictions. Therefore, it appears difficult 
for dynamo theories to produce the fields observed 
on TTSs. As a result, we suggest that the fields 
on these young stars may indeed be primordial in 
origin, as suggested by Tayler (1987) and Moss 
(2003). 
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Table 1 
Sample of Stars 





Spec. 


M* 




M x 10 s 


Prot 


Beq 


B e q 


B eq 


B eq 


L x 


Star 


Type 


(M ) 


(R©) 


(Mgyr- 1 ) 


(days) 


(G) 


(G) 


(G) 


(G) 


(10 30 erg s" 1 ) 


AA Tau 


K7 


0.70 


1.74 


0.33 


8.20 


1020 


290 


1210 


380 


0.439 


BP Tau 


K7 


0.70 


1.99 


2.88 


7.60 


1850 


620 


2180 


880 


0.707 


CY Tau 


Ml 


0.40 


1.63 


0.75 


7.90 


1130 


380 


1330 


530 


< 0.627 


DE Tau 


M2 


0.32 


2.45 


2.64 


7.60 


490 


190 


580 


230 


< 0.287 


DF Tau 


Ml 


0.40 


3.37 


17.7 


8.50 


680 


290 


790 


320 


0.646 


DG Tau 


K7.5 


0.65 


2.05 


4.57 


6.30 


1610 


560 


1900 


760 


< 0.223 


DH Tau 


M1.5 


0.36 


1.39 


0.10 


7.20 


550 


160 


640 


260 


1.318 


DK Tau 


K7 


0.68 


2.49 


3.79 


8.40 


1190 


410 


1400 


560 


< 0.223 


DN Tau 


MO 


0.51 


2.09 


0.35 


6.00 


320 


100 


380 


150 


0.564 


GG Tau A 


K7 


0.68 


2.31 


1.75 


10.3 


1280 


420 


1510 


610 


< 0.107 


GI Tau 


K6 


0.93 


1.74 


0.96 


7.20 


1900 


560 


2240 


900 


0.241 


GK Tau 


K7 


0.69 


2.15 


0.64 


4.65 


390 


110 


450 


180 


0.241 


GM Aur 


K7 


0.70 


1.78 


0.96 


12.0 


2540 


800 


2990 


1200 


0.562 


T Tau 


KO 


2.30 


3.31 


4.40 


2.80 


420 


110 


490 


200 


0.568 


TW Hya 


K7 


0.74 


1.00 


0.20 


2.20 


950 


250 


1120 


450 


1.380 



Note. — All predicted fields are the equatorial magnetic field strength for an assumed dipolar magnetic field. In 
order, the predicted fields come from Konigl (1991), Collier Cameron and Campbell (1993), Shu et al. (1994), and 
Long ct al. (2005). 
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Table 2 
Observing Runs 





Obs. 


Wavl. 


Star 


Run 


Set 


AA Tau 


III 


1,2,3 


BP Tau 


II 


1,2,3 


CY Tau 


IV 


1,2,3 


DE Tau 


III 


1,2,3 


DF Tau 


II 


1,2,3 


DG Tau 


III 


1,2,3 


DH Tau 


III 


1,2,3 


DK Tau 


II 


1,2,3 


DN Tau 


III 


1,2,3 


GG Tau A 


III 


1,2,3 


GI Tau 


III 


1,2,3 


GK Tau 


III 


1,2,3 


GM Aur 


I 


2 


T Tau 


II 


1,2,3 
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Table 3 
Results 







B 


B eq 


Lx 


Star 




(kG) 


(kG) 


(10 30 erg s- 1 ) 


AA Tau 


0.52 ±0.09 


2.78 


1.02 


4.32 


BP Tau 


1.08 ±0.25 


2.17 


1.03 


4.42 


CY Tau 


0.65 ±0.21 


1.16 


1.04 


1.36 


DE Tau 


1.16 ±0.17 


1.12 


1.04 


3.34 


DF Tau 


1.37 ±0.15 


2.90 


1.04 


20.7 


DG Tau 


2.67 ±0.41 


2.55 


1.00 


8.19 


DH Tau 


1.77 ±0.49 


2.68 


1.42 


2.20 


DK Tau 


1.46 ±0.23 


2.64 


1.02 


9.27 


DN Tau 


0.54 ±0.09 


2.00 


1.01 


4.51 


GG Tau A 


0.35 ±0.19 


1.24 


1.02 


3.28 


GI Tau 


1.00 ±0.11 


2.73 


1.44 


4.23 


GK Tau 


1.47 ±0.31 


2.28 


1.02 


5.59 


GM Aur 


1.25 


2.22 


1.02 


3.51 


T Tau 


0.45 ±0.77 


2.37 


0.84 


15.8 


TW Hya a 


0.07 ±0.04 


2.61 


1.76 


1.13 



a Veiling and field measurements taken from Yang ct al. (2005). 
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Table 4 
Range of Predicted Fields 





VBJ 




HEG 




GHBC 




CG 




WG 






M x 10 8 


Be," 


M x 10 8 


Be," 


M x 10 8 


-Beg* 1 


M x 10 8 


Beg" 1 


M x 10 8 




Star 


(Moyr- 1 ) 


(kG) 


(Moyr" 1 ) 


(kG) 


(Moyr" 1 ) 


(kG) 


(Moyr" 1 ) 


(kG) 


(Moyr- 1 ) 


(kG) 


AA Tau 


0.71 


1.56 


12.6 


4.04 


0.33 


0.96 


0.40 


1.05 


0.65 


2.34 


BP Tau 


2.43 


2.46 


15.8 


4.06 


2.88 


1.62 


2.30 


1.45 


1.32 


1.81 


CY Tau 






0.63 


2.62 


0.75 


1.38 


0.80 


1.43 


0.14 


0.75 


DE Tau 


18.0 


0.07 


31.6 


1.18 


2.64 


0.49 


4.60 


0.65 


4.07 


0.62 


DF Tau 


21.9 


0.37 


125.9 


0.67 


17.7 


0.57 


20.8 


0.62 


1.00 


2.03 


DG Tau 






199.5 


9.09 






91.9 


7.36 


4.57 


2.44 


DH Tau 


2.83 


0.25 














0.11 


0.93 


DK Tau 


0.61 


1.93 


39.8 


2.19 


3.79 


0.95 


0.50 


0.35 






DN Tau 


0.13 


0.15 


3.20 


0.84 


0.35 


0.30 


0.20 


0.23 


0.19 


0.38 


GG Tau A 


3.01 


1.48 


20.0 


1.41 


1.75 


1.05 


2.50 


1.26 


1.26 


1.79 


GI Tau 


2.02 


9.16 


12.6 


1.06 


0.96 


1.70 


0.60 


1.35 


0.83 


1.19 


GK Tau 


0.05 


0.23 


6.30 


0.86 


0.64 


0.32 


1.40 


0.48 


0.65 


0.57 


GM Aur 


0.74 


4.30 


2.50 


5.20 


0.96 


2.34 


0.70 


2.00 


0.66 


6.86 


T Tau 


















3.16 


0.39 


Cor b 


0.27 


0.43 


0.21 


0.51 


0.07 


0.85 


0.21 


0.52 


0.18 


0.55 



a The predicted field strength is the equatorial field strength of an assumed dipolc field using the Shu ct al. (1994) model. 

b Given on this line in each M column is the value of the correlation coefficient, r, between the predicted field strenghts and the observed mean 
field strengths. Given on this line in each B eq column is the associated false alarm probability. 
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Table 5 
Correlations of Derived 
Magnetic Parameters with 
Stellar Parameters 



Quantities 



Compared 


r 


f P 


B vs P- Q \ 


0.22 


0.44 


$ vs P ro ] 


0.05 


0.85 


B vs t c P- \ 


0.23 


0.41 


$ VS T c P rQ \ 


-0.24 


0.40 


B vs t c 


0.03 


0.92 


$ vs r c 


-0.57 


0.03 


B vs T eff 


0.26 


0.36 


$ vs T eff 


0.37 


0.18 


B vs M* 


0.42 


0.12 


$ vs 


0.20 


0.48 


i? vs L* 


0.10 


0.73 


$ vs L* 


0.66 


< 0.01 


i? vs Age 


0.23 


0.41 


$ vs Age 


-0.31 


0.26 



Note. — Convective turnover times, 
t c , were kindly computed by Y.-C. Kim, 
based on models described in Kim & 
Demarque (1996) and Preibisch et al. 
(2005). 
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Fig. 1. — K band spectra of DK Tau are shown 
in the black histogram. The top 2 panels show 
Zeeman sensitive Ti I lines and a Sc I line. The 
bottom panel shows Zeeman insensitive CO lines. 
The dashed curve shows a model with no magnetic 
field. The smooth solid curve shows the model fit 
for a star covered by a mixture of regions which 
include 0, 2, 4, and 6 kG fields. Here, the mean 
field strength averaged over the stellar surface is 
2.64 kG. 
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1.05 




2.220 2.221 2.222 2.223 2.224 

Wavelength (|im) 




P.as H , , 

2.227 2.228 2.229 2.230 2.231 

Wavelength (|im) 



1.05 




pus h ..... 

2.310 2.311 2.312 2.313 2.314 
Wavelength (|im) 



Fig. 2. — K band spectra of DE Tau are shown 
in the black histogram. The top 2 panels show 
Zeeman sensitive Ti I lines and a Sc I line. The 
bottom panel shows Zeeman insensitive CO lines. 
The dashed curve shows a model with no magnetic 
field. The smooth solid curve shows the model fit 
for a star covered by a mixture of regions which 
include 0, 2, 4, and 6 kG fields. Here, the mean 
field strength averaged over the stellar surface is 
1.12 kG. 
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Fig. 3. — The observed mean magnetic field 
strength, B, plotted versus the predicted equato- 
rial magnetic field strength from the Shu et al. 
(1994) treatment of the magnetospheric accretion 
model. 



Fig. 4. — The measured X-ray luminosity plot- 
ted against the predicted X-ray luminosity for our 
sample of stars. The line of equality is shown by 
the dashed line. The majority of these CTTSs 
appear underluminous in X-rays given their mag- 
netic field properties. 
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Observed Mean Magnetic Field (kG) 




Fig. 5. — The observed mean magnetic field on 
our sample of CTTSs plotted versus the maximum 
allowed magnetic field strength predicted assum- 
ing pressure equilibrium with a surrounding non- 
magnetic atmosphere. The observed fields tend to 
be much stronger than the pressure equilibrium 
predictions, implying magnetic filling factors near 
unity. 
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